Kink-induced transport and segregation in oscillated granular layers 
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We use experiments and molecular dynamics simulations of vertically oscillated granular layers to 
study horizontal particle segregation induced by a kink (a boundary between domains oscillating out 
of phase) . Counter-rotating convection rolls carry the larger particles in a bidisperse layer along the 
granular surface to a kink, where they become trapped. The convection originates from avalanches 
that occur inside the layer, along the interface between solidified and fluidized grains. The position 
of a kink can be controlled by modulation of the container frequency, making possible systematic 
harvesting of the larger particles. 

PACS numbers: 45.70.Mg, 81.05.Rm, 44.27.+g, 05.60.-k 



In 1831 Faraday observed convective motion of grains 
in heaps in vertically oscillated granular layers jlj . Later 
there weie many studies of convection in oscillated gran- 
ular materials |a,llllll0,Il,li,llESElE3 These 
convection phenomena are driven by interstitial air |^ 0| 
and shear due to sidewalls [7], llOl lll| . One consequence 
of the sidewall-driven convection is vertical segregation 
of grains of different sizes, which is an example of the 
Brazil-nut effect [HEIEI- 

We consider here a different kind of convection, one 
that is driven not by interstitial air or sidewalls but rather 
arises from the intrinsic dynamics of the oscillated gran- 
ular layers. This convection is associated with kinks, 
which are boundaries separating oscillating domains of 
opposite phase; when the layer on one side of a kink is 
moving up, the layer on the other side of the kink is mov- 
ing down. Kinks can form in oscillating layers when the 
maximum container acceleration a^ax is large enough 
[o-max > 4.5^, where g is the gravitational acceleration) 
so that the granular layer hits the container bottom only 
every other cycle |l6l K^. Kinks spontaneously form in 
a layer for Umax > 7g U^- 

During each oscillation cycle, the granular material on 
one side of a kink strikes the container bottom and solidi- 
fies, while the fluidized grains on the other side avalanche 
down along the front of the solidified grains. We exam- 
ine this avalanche process, which will be shown to be 
responsible for convection associated with kinks. In a 
bidisperse layer, such convective motion conveys larger 
grains toward the kink. We will show that controlled 
motion of a kink leads to horizontal size segregation, in 
contrast to the vertical size segregation in the Brazil-nut 
effect. We will first describe the laboratory observations, 
and then show how the particle transport, trapping, and 
horizontal segregation can be understood using molecular 
dynamics simulations. 

Experimental observations: Trapping, transport, and 
segregation — We find that larger grains in the vicinity 
of a kink in an oscillating bidisperse layer move toward 
the kink and remain there, while the larger grains far 
from a kink move only diffusively, as the solid white lines 




FIG. 1: Segregation and controlled transport of 650 fim di- 
ameter black glass spheres by a kink in an experiment with a 
10-particle deep oscillated layer of 165 /xm bronze spheres at 
Omax ~ 5g and / = 92 Hz. The motion of the kink is con- 
trolled by modulation of the phase difference A(j> between the 
primary oscillation signal and an added small subharmonic 
sinusoidal perturbation (see text), (a) A kink (dashed line) 
sweeps across the layer after a rapid change in A<j) of 2n/3; 
glass spheres close to the kink move toward it and remain 
trapped in the kink at the surface of the layer, while far from 
the kink they diffuse; these two types of trajectories are illus- 
trated by the solid white lines. The numbers in each frame 
denote plate oscillations after A(j} is changed, (b) Sinusoidal 
oscillation of a kink, with A(j}{t) = fms/fmrSm{2Tifmrt), 
fms = 0.17 Hz and fmr =0.1 Hz. The large trapped par- 
ticles follow the motion of the kink. 
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FIG. 2: Grains form a pair of convection rolls, flowing down- 
ward at the kink, as illustrated by these projections of a layer 
with a kink at different times ft during a cycle {ft — when 
the container is at its equilibrium position, moving upward), 
obtained from a simulation {amax ~ 5.2g, f — 69 Hz). The 
vertical position of the center of mass as a function of time of 
the left (right) side of the kink is shown as dashed (solid) line 
in the inset of the top panel. The two bold vertical arrows 
on either side of the kink indicate the directions of motion of 
these two domains. The small arrows show the discrete grain 
velocity u{r{t),t) averaged over the shorter horizontal direc- 
tion. Each circle represents a grain, and the container bottom 
is indicated by horizontal gray lines. The region shown is far 
from the rigid sidewalls in the longer horizontal direction. 



in Fig. n (a) illustrate. Further, if the location of a kink 
is controlled by modulating the container oscillation fre- 
quency |0| , the large grains follow the motion of the kink 
(Fig. n (b)); this leads to a horizontal size segregation. 
We consider a 10-particle deep layer of 165 /j,m bronze 
spheres (mass density 8.3 g/cm^) together with a few 
hundred of 650 /im glass spheres (2.5 g/cm"^), in an evac- 
uated square container of horizontal area 8.9 x 8.9 cm^ 
(the pressure is 5 Pa). The layer is subject to a verti- 
cal sinusoidal oscillation with an amplitude A and a fre- 
quency /; the oscillation is characterized by two of the 
following three parameters: the maximum acceleration 
amax = ^(27r/)^, the maximum velocity Vmax = 27rA/, 
and the frequency /. In some experiments we add a 
perturbative subharmonic forcing to control the posi- 
tion of a kink, where the plate position z is given by 
z{t) = Asin(27r/t) + 0.01Asin(7r/i -I- A(t>{t)). This sec- 
ondary forcing breaks the symmetry of the momentum 
transfer across the kink, which makes the kink move [l8| . 



The relative phase between the main oscillation and the 
subharmonic perturbation controls the position of the 
kink within the container (see Fig. 

We consider the range 4.5g < Gmax < 5.5g, where the 
layer remains flat except for the kinks; similar transport 
and segregation occur for larger amax where wave pat- 
terns form, but it is simpler to quantify the transport 
and segregation processes for a non-patterned flat layer. 

Microscopic dynamics; kink-associated convection — 
The experimental observations reveal that a mean flow 
toward a kink exists in the vicinity of the kink. We use 
a previously validated three-dimensional (3-d) molecular 
dynamics simulation 0, to understand the trap- 
ping and transport produced by the microscopic dynam- 
ics. The coefficient of restitution e and the coefficient 
of friction /x were set to 0.7 and 0.5 respectively, as in 
Ref. [i^l. We simulate 8-particle deep monodisperse lay- 
ers of spherical grains of diameter 165 fim in a rectangular 
box of horizontal area 160cr x lOcr, where a is the particle 
diameter. Periodic boundary conditions are used in the 
shorter horizontal direction (y-axis), and vertical rigid 
sidewalls are imposed in the longer horizontal direction 
(x-axis). The rigid sidewalls and the bottom plate are 
assumed to be made of the same material as grains. 

We distinguish two types of velocities, the instanta- 
neous grain velocity v(r(t),t) = fiinAt^o+(r(i + Ai) — 
r{t))/At and the discrete grain velocity u(r(t),i) — 
{r{t + 2T) - r(t))/(2r), where r{t) is the position vector 
of a grain at a given time t, and 2T = 2// is the oscilla- 
tion period of the macroscopic state. The discrete veloc- 
ity corresponds to the grain velocity in the //2-strobed 
frame. Following the trajectories of individual grains, we 
observe that grains do not come back to the same posi- 
tions after 2T, even though the layer returns to the same 
macroscopic state. Instead, in the discrete velocity field, 
grains form a pair of counter-rotating convection rolls 
that flow downward at the kink (Fig. ^ . The shape of 
the rolls changes during a cycle, as does the shape of the 
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FIG. 3: A schematic diagram showing trapping of larger 
particles at a kink in a layer of smaller particles. The larger 
particles move toward the kink due to the surface flow of 
the convective motion of the smaller particles. The larger 
particles do not flow down in between the convection rolls, 
but remain there, due to their size. Dashed arrows represent 
the convection rolls formed by the smaller particles, and solid 
arrows indicate the direction of the motion of large particles. 
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FIG. 4: The trapping zone width (defined in the inset) for 
various amax's and /'s, obtained in simulations of 8-particle 
deep monodisperse layers of 165 fim spheres. The width is 
linearly proportional to Vmax, as shown by the dashed line (a 
least square fit). The symbols correspond to: •, Umax ~ 4.9(7; 
o, 5g; *, 5.1g; □, 5.2g; A, 5.3g. For each Umax, the data 
points correspond to / = 86 Hz, 78 Hz, 69 Hz, 60 Hz, 52 
Hz, respectively. Inset: The trapping zone is defined as the 
region where the vertically integrated |up (gray solid line, in 
an arbitrary unit) is larger than some small value (dash-dotted 
line) . 



layer, but the sense of circulation remains the same. For 
the parameters of Fig.|2 it takes about lOOT for grains in 
the outer- most part of the rolls to complete a full circuit. 

Trapping and transport — The presence of the convec- 
tion rolls explains the motion of large grains toward a 
kink: Large grains far from the kink rise to the surface 
as a consequence of the Brazil-nut effect and those that 
enter the convecting region will move toward the kink 
due to the surface flow of the convective motion (Fig.|21l. 
When the large grain reaches the down-flow channel re- 
gion between the convection rolls, it will not flow down 
because it is too large, and it becomes trapped there, as 
illustrated in Fig. ^ We have observed trapping for a 
wide range of impurity particle diameters and densities 
relative to those of the oscillated bronze spheres. We 
introduced hundreds of the same type of impurity parti- 
cles, including lead, stainless steel, glass, and polystyrene 
spheres of diameter ratios 1.12 to 11.8, density ratios 0.14 
to 1.4, and mass ratios 0.4 to 2500. Trapping occurred 
in all cases studied except for the diameter ratio below 
1.12, where the motion of the impurity grains was vi- 
sually indistinguishable from the motion of the bronze 
spheres. 

We refer to a region in which large grains move toward 
a kink and become trapped as the "trapping zone" . The 
width of the trapping zone, defined in Fig. 01 is found to 
be linearly proportional to Vmax for the range of oscilla- 
tion parameters in our study (Fig. ^J. There is a limit 
to the number of large grains that can be trapped; for 




FIG. 5: Simulations reveal an avalanche occurring inside 
the layer, along the interface between the solidified and fiu- 
idized parts: (a) Instantaneous grain velocities v{r{t),t) in 
the laboratory frame (indicated by arrows), averaged over the 
shorter horizontal direction, when the left side of a kink is be- 
ing pushed up by the bottom plate (a horizontal solid line) 
for the same case as in Fig. |21 A gray background repre- 
sents the layer, (b) Instantaneous velocities in the container 
frame, during the same time. As the left side is pushed up, 
it solidifies on the container bottom, while the right side is 
fluidized and cascades down onto the left side, as shown by 
gray arrows, {umax ~ 5.2g, f — 69 Hz). 



the range of Vmax studied in this paper, the maximum 
number of trapped particles corresponds to the number 
of large grains forming a monolayer above the trapping 
zone. In the case of Fig. ^ (a), Vmax of the primary os- 
cillation is 0.087 m/s, and the trapping zone width from 
the last frame (t = 566T) is seen to be about 5 glass 
particle diameters or 20 bronze particle diameters, which 
is comparable to the simulation result in Fig. 01 

As Vmax increases, qualitative changes occur in the 
motion of large grains: The down- flow channel region di- 
lates more (i.e., the granular volume fraction decreases) 
so that there is enough room for large grains to circulate 
with the convecting smaller grains throughout the vol- 
ume of the trapping region, disappearing in the kink and 
popping up from the upward flow region of the convec- 
tion rolls. The maximum number of trapped particles in 
this regime is not determined by the area of the trapping 
zone; we have not investigated this regime. 

We observe in the experiment that when a kink moves 
too rapidly, large grains leak out of the kink; there is a 
maximum speed of a kink for which it can convey large 
grains. The mechanism in Fig. |31 suggests that if a kink 
moves more than the trapping zone width during a cy- 
cle, the large grains would not follow the motion; the ex- 
perimental observations are consistent in that when the 
speed of a kink exceeds approximately twice the layer 



4 



depth (comparable to the zone width in Fig. ^ per time 
2T, large grains are left behind the kink. 

Internal avalanches and convection — An examination 
of the instantaneous velocity field v{r{t),t) reveals that 
the convection is driven by avalanches oscillating per- 
pendicular to the kink. When the left side of a kink is 
pushed up, the right side cascades down onto the left side 
of the kink (Fig.EKa)). These avalanches inside the layer 
are visualized better when the velocity field is viewed in 
the container frame (Fig. (b)): When the left side is 
pushed up, it solidifies and moves with the same veloc- 
ity as the container until it takes off from the container. 
Meanwhile, the right side is still falling and is fluidized. 
This leads to a density gradient between the solidified 
and the fluidized regions. Grains flow along the surface 
of the solidified region, toward the right side. We refer to 
this motion inside the layer as an internal avalanche, to 
contrast with the usual avalanche, which occurs on a sur- 
face • (Such an internal avalanche was suggested by 
Laroche et al. |^ as a mechanism for convection in heap- 
ing.) In the following cycle, another internal avalanche 
occurs in the opposite direction, pushing grains back to 
the left side and to lower heights. The oscillatory inter- 
nal avalanches result in the downward flow at the kink. 
At impact with the container, grains at the bottom of 
the layer move away from the kink in a horizontal direc- 
tion and later rise, forming the upward flow at the outer 
side of the rolls. The trapping zone width (or convection 
roll size) is determined by the horizontal displacements 
of such grains. 

The convective motion that has been observed previ- 
ously in subharmonic standing wave patterns in oscil- 
lated layers jl^ can be explained by the same mechanism: 
when a patterned layer collides with the container, large 
density gradients form along the pattern inside the layer. 
Internal avalanches perpendicular to these gradients lead 
to convection rolls. 

Conclusions — We have described a form of convection 
in granular media that is driven not by interaction of the 
grains with air or the container sidewalls but by the in- 
trinsic dynamics of the layer when a kink is present. The 
kink and internal avalanches that drive this convective 
motion are unique to oscillating granular media because 
the solidification and fluidization during each cycle that 
occur in the granular layer do not occur in an oscillated 
liquid layer. 

The kinks that lead to the avalanche motion and con- 
vection are ubiquitous in oscillated granular layers at 
high enough container accelerations {omax > ^g); hence 
kink-associated convection and trapping of large parti- 
cles are generic features of bidisperse oscillated granular 
layers in that regime. Since the location of a kink can 
be controlled by modulation of the container oscillation 



frequency, it is possible to harvest the segregated larger 
particles in a controlled way by sweeping a kink to the 
edge of a container. 

This work was supported by the Engineering Research 
Program of the Office of Basic Energy Sciences of the 
U. S. Department of Energy (Grant No. DE-FG03- 
93ER14312) and The Texas Advanced Research Program 
(Grant No. ARP-055-2001). 



* "moon@chaos.utexas.edu' 
t swift@chaos.utexas.edu 
[1] M. Faraday, Phil. Trans. R. See. London 52, 299 (1831). 
[2] W. Kroll, Chem. Ingr. Tech. 27, 33 (1955). 
[3] G. Ratkai, Powder Technol. 15, 187 (1976). 
[4] S. B. Savage, J. Fluid Mech. 194, 457 (1988). 
[5] C. Laroche, S. Douady, and S. Fauve, J. Phys. France 

50, 699 (1989). 
[6] H. K. Pak, E. Van Doom, and R. P. Behringer, Phys. 

Rev. Lett. 74, 4643 (1995). 
[7] E. Clement, J. Duran, and J. Rajchenbach, Phys. Rev. 

Lett. 69, 1189 (1992). 
[8] J. A. C. Gallas, H. J. Herrmann, and S. Sokolowski, Phys. 

Rev. Lett. 69, 1371 (1992). 
[9] Y-h. Taguchi, Phys. Rev. Lett. 69, 1367 (1992). 
[10] J. Lee, J. Phys. A 27, L257 (1994). 

[11] J. B. Knight, H. M. Jaeger, and S. R. Nagel, Phys. Rev. 
Lett. 70, 3728 (1993); E. E. Ehrichs, H. M. Jaeger, G. 
S. Karczmar, J. B. Knight, V. Yu. Kuperman, and S. 
R. Nagel, Science 267, 1632 (1995); J. B. Knight, Phys. 
Rev. E 55, 6016 (1997); E. L. Grossman, ibid. 56, 3290 
(1997). 

[12] R. D. Wildman, J. M. Huntley, J.-P. Hansen, D. J. 

Parker, and D. A. Allen, Phys. Rev. E 62, 3826 (2000); 

R. D. Wildman, J. M. Huntley, and D. J. Parker, Phys. 

Rev. Lett. 86, 3304 (2001). 
[13] R. L. Brown, J. Inst. Fuel 13, 15 (1939). 
[14] J. C. Williams, Powder Technol. 15, 245 (1976). 
[15] M. E. Mobius, B. E. Lauderdale, S. R. Nagel, and H. M 

Jaeger, Nature, 414, 270 (2001). 
[16] F. Melo, P. B. Umbanhowar, and H. L. Swinney, Phys. 

Rev. Lett. 75, 3838 (1995). 
[17] P. B. Umbanhowar, F. Melo, and H. L. Swinney, Physica 

A 249, 1 (1998). 
[18] S. J. Moon, M. D. Shattuck, C. Bizon, D. I. Goldman, 

J. B. Swift, and H. L. Swinney, Phys. Rev. E 65 011301 

(2002). 

[19] I. S. Aranson, D. Blah, W. K. Kwok, G. Karapetrov, U. 

Welp, G. W. Crabtree, V. M. Vinokur, and L. S. Tsim- 

ring, Phys. Rev. Lett. 82, 731 (1999). 
[20] C. Bizon, M. D. Shattuck, J. B. Swift, W. D. McCormick, 

and H. L. Swinney, Phys. Rev. Lett. 80, 57 (1998). 
[21] J. Rajchenbach, Phys. Rev. Lett. 89, 074301 (2002). 
[22] C. Bizon, M. D. Shattuck, J. R. de Bruyn, J. B. Swift, 

W. D. McCormick, and H. L. Swinney, J. Stat. Phys. 93, 

449 (1998). 



